Triplet pairing due to spin-orbit-assisted electron-phonon coupling. 
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We propose a microscopic mechanism for triplet pairing due to spin-orbit- assisted electron inter- 
action with optical phonons in a crystal with a complex unit cell. Using two examples of electrons 
with symmetric Fermi surfaces in crystals with either a cubic or a layered square lattice, we show 
that spin-orbit-assisted electron-phonon coupling can, indeed, generate triplet pairing and that, in 
each case, it predetermines the tensor structure of a p-wave order parameter. 

PACS numbers: 74.20.Fg, 74.20.Rp, 74.20.-z, 74.25.Kc, 74.70.Pq 
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Spin-triplet superconductivity has been found in sev- 
eral materials including SY2-RWO4 0] and uranium-based 
compounds (UPt 3 , UBe 13 , UGe 2 , URuGe) HHH@. 
While the phenomenology of the triplet state is well 
known (it is similar to that of p-wave-superfluid He 3 0, 
7]), the microscopic origin of the triplet pairing in each of 
these materials and its relation to the band structure are 
the subject of intense investteationjl HH 1100, EH 03, 
and the existing models [50 Hllll2l | attribute the nature 
of the electron-electron attraction in the triplet channel 
to the interaction via magnetic fluctuations. 

In this Letter, we show that a strong spin-orbit cou- 
pling in constituent atoms of a metal with a complex 
unit cell can generate triplet pairing via mechanism 
which involves optical phonons (instead of spin waves). 
The attention 13] to spin-orbit coupling in the elec- 
tron band structure in the context of odd-parity pair- 
ing has been attracted by the observation of supercon- 
ductivity in the inversion-asymmetric CePt 3 Si |l4|. In 
other known triplet superconductors, spin-orbit coupling 
in the band structure can be ignored due to the inver- 
sion symmetry of their crystalline lattice. However, in 
center-symmetric crystals with a multi-component unit 
cell sublattice displacements may break inversion symme- 
try, so that atomic spin-orbit coupling can give rise to an 
additional, spin-dependent interaction between electrons 
and optical phonons. Such a spin-orbit-assisted electron- 
phonon (e-ph) coupling arises from a spin-dependent en- 
velope between orbitals belonging to atoms from different 
sublattices (mutially displaced by the vibration in the op- 
tical phonon mode), and it may generate attraction in the 
spin-triplet Cooper channel via the BCS-likc mechanism. 

Below, we demonstrate that the spin-orbit-assisted e- 
ph coupling does, indeed, lead to spin-triplet p-wave pair- 
ing. The symmetry of the order parameter of the result- 
ing superconducting state is determined by the particular 
crystalline lattice and the polarization properties of the 
strongest-coupled optical phonon modes, which we show 
for two model examples: of (i) a cubic lattice with two 
atoms in the unit cell and a spherical Fermi surface and 
(ii) a layered structure with in-plane square symmetry 
and a cylindrical Fermi surface. 



For a cubic center-symmetric two-component lattice 
the electron spin coupling to the inversion symmetry 
breaking sublattice displacement w can be described phe- 
nomenologically using the invariant [l5| 



V c = (a/2) (w- [p x a] + [p x a] -w) . 



(1) 



In the second-quantized form, this can be represented as 



(2) 



V=J2 ^i(P»«j)^p+f/«p-} («qC+«- q< 
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V c c = (a/ v / 2uj c m c ) l c [p x &] , 

Here, a) (a) and f'(f) stand for creation(annihilation) 
operators of phonons and electrons, 1^, to^ and are 
the polarisation vector, frequency and reduced mass of 
the optical phonon £, indices fi,v identify the electron 
in the Kramers doublet for each momentum (spin for a 
simple conduction band), and the parameter a is a ma- 
terial constant. In a cubic crystal, we neglect a small 
difference between the LO and TO phonon frequencies, 
thus k w . In a less symmetric crystal, V** = 

T,k=x, y ,zV kA k(p><i)> where ^fc(p.-q) = 4(p>q) and 

4(-P,q) = -AM, t^s, tK(-p,q) = -?<(p,q), 
V (Pi — l) = » (Pi q)- F° r example, in a layered crystal 
with a square lattice where we neglect the electron mo- 
tion perpendicular to the x-y planes, we take into account 
two invariants analogous to Eq. |Q, a±w z (p x a v — p y a x ) 
and aytj 2 (w x p y — w v p x ), to distinguish between the elec- 
tron coupling to optical phonons polarized perpendic- 
ular (C = z;uj z = lo±) and in the plane of the layers 

(C = %,y;v x ,y = w||): 
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A spin-dependent e-ph coupling, Eq. El and a virtual 
polarisation of the crystalline unit cell |lfij ] result in the 
effective spin-dependent interaction between electrons on 
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the opposite sides of the Fermi surface, 
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(4) 
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In Eq. (JSJ, we neglect the dispersion of optical phonon 
modes, w^(q) ~ u>£ and take into account the quasipar- 
ticle states close to the Fermi surface (p = |p| » pf), 
with |e Pi | < cj^, where e p is the quasiparticle energy 
calculated with respect to the Fermi level. The effec- 
tive electron-electron (e-e) interaction, U is spin- and 
momentum-dependent, therefore, it may allow for attrac- 
tion in the spin-triplet channel and odd pairing. 

To determine what triplet phase can be formed in 
a metal with the Cooper interaction in Eas. (|4l5fl . one 
needs to identify particle-particle channels in which the 
electron-electron interaction is attractive. The classifica- 
tion of pairing channels and the resulting triplet order 
parameter is usually based upon the symmetry of the 
electronic Hamiltonian Q . Following a procedure devel- 
oped for p-w&ve superfluidity in He 3 0, , we determine 
operators of singlet (b) and triplet (t s ) electron pairs, 

Hp) = irf v U-pfvp; 4(p) = i{^ v ^ s )^f^ P fu P - (6) 

In terms of these operators 0|i the Hamiltonian U in 
Eq. Q can be written down in a separable form, 

U = E [^ SlS2 4 1 (Pi)^(P2) + CSt(p 1 )6(p 2 )] , (7) 
P1P2 

which includes the interaction both in the singlet chan- 
nel [lg, with the constant C = M"$ 'ci^er^, and in the 
triplet channel, with 
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Below, we limit the analysis to the case of p-wave pair- 
ing (excluding higher-order angular harmonics) described 
by the tensor order parameter, 



B ms = y^n m (4(p)); 
p 



n = — 



(9) 



and focus it on two model examples: (i) a 3D metal 
with a cubic lattice and a spherical Fermi surface and 
(ii) a layered material with a square in-plane lattice and 
a cylindrical Fermi surface. In such symmetric systems, 
triplet pairs can be characterised by well-defined quan- 
tum number related to irreducible tensor representations 
of the corresponding symmetry group. 

(i). In the case of a model metal with a cubic lattice 
and a spherical Fermi surface the symmetry is the full 
rotational group SO3. The involment of the spin-orbit 



coupling in the formation of the pair interaction assumes 
that the spin and orbital degrees of freedom of the Cooper 
pair are not independent, so that its total spin 5 = 1 and 
orbital angular momentum, L = 1 sum into the total 
angular momentum, J = 0, 1, 2. The tensor structure of 
the order parameter for these values of J is the following: 
B^=° = 8 ms De^, where B is a scalar; B^ 1 = e msl D h 
where D is a complex vector; and B^ 2 is a traceless 
symmetric tensor. 

Thus, decompose e-e interaciton in Eq. into three 
channels corresponding to different angular momenta of 
the pair. Using the relation of the matrix M"j? with the 

e-ph coupling V£ in Eq. ©, we find that 
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The first term in the matrix G^J, 2 is a product of pro- 
jectors onto the J = paired state; the second projects 
onto the states with J = 1, whereas the last term selects 
only the states with J = 2. Interaction in each channel 
is described by the corresponding constant, 



9J = xjg so , g so ~ (a/ui ) 2 p 2 F /2m 



(10) 



where the scale of the interaction depends on the strength 
of the original electron-phonon interaction for electrons 
near the Fermi surface, frequency lu and reduced mass 
to of the optical phonon (here we disregard a small 
difference between longitudinal and transverse optical 
phonon modes in a cubic crystal). Interaction in chan- 
nels with different J's differs by the numerical factors, 



(11) 



According to Eq. JJJ), the positive sign of gj corre- 
sponds to attraction. Since only one of the effective cou- 
pling constants gj — g so xj is positive, we conclude that 
in the cubic system Cooper interaction due to the spin- 
orbit-assisted e-ph coupling leads to the condensate of 
pairs with the total angular momentum J = 1. 

By selecting a pairing channel with J — 1 we do 
not fully determine the order parameter B ms in the 
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lowest energy triplet state, yet. In B£ 

D =De u ? (1 + il') /\J\ 2 + (l') 2 is a complex vector 
parametrized using two perpendicular vectors, 1 and 1', 
l 2 > (!') . To find a favourable state, we minimize the 
condensate energy, 



^ = E(I £ pI- £; pm)+E^I s " 



(12) 
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TABLE I: J = 1 phases due to spin-orbit-assisted electron- 
phonon coupling in a metal with spherical Fermi surface and 
cubic lattice. 



Vector D 


mean-field D 


Energy, £ 


le^D 


D = {2L0 o /g)e-^ 9V 


£ = -vg 2 D 2 /3 




D = {V2uo /g)e- & ' 9V 


£ = -vg 2 D 2 /G 



as a function of D, which includes finding the relative 
size of two vectors 1 and 1' and the absolute value of D. 
The expression for £ in Eq. H12[l can be formally ob- 
tained using the Hubbard-Stratonovich transformation 
which splits the electron-electron interaction in Eq. fus- 
ing the held B ms and, then, by integrating out fermionic 
degrees of freedom. Equation is the self-consistency 
equation for the minimum of £ , and 



|d P | 2 



[dp 



(13) 



describes the spectrum of excitations 0, in the phase 
with order parameter B ms , where three components of 



the complex vector d p are defined as 



,B 



m-'-'ms • 



(14) 



The vector product id p x d* determines the spin quanti- 
zation axis and spin splitting for Bogolyubov quasiparti- 
cles for each direction n = p/p on the Fermi sphere and 
fi = ±1 denotes the polarisation of the quasiparticle spin 
with respect to such an axis. In Eqs. (|12I13|) . e p is the 
quasiparticle energy in the normal state calculated with 
respect to the Fermi level, and g is the interaction con- 
stant in the attraction channel (here, g — g\ = >cig so ). 

In the context of pairing with J = 1 and B^ 1 = 
e msl Di, two possible extrema of the function £ can be 
found using the analysis of the degeneracy space of the 
order parameter |lj|. For fixed |D| = D, £(D) takes 
extremal values on vectors D that cover minimal sur- 
faces under the application of SO3 rotations. Two such 
possibilities arc Do= le^D where 1 is a real unit vector 
(that's, 1' = 0), and Di=^e^£>, where 1 and 1' are 
two perpendicular unit vectors. For each of these two 
cases we find the value of D by minimizing the ground 
state energy £ in Eq. i|12|) . The results |20J are shown 
in Table ^indicating that the phase with Do= le^D and 
an anisotropic gap vanishing along the direction 1 (since 
d = gD[n x 1]) is energetically favorable. 

Although the example of a cubic crystal with highly 
symmetric Fermi surface may be too abstract, neverthe- 
less, it illustrates that the spin-orbit-assisted e-ph cou- 
pling can generate triplet pairing between electrons and 
that the suggested pairing mechanism is highly prescrip- 
tive for the structure of the resulting order parameter. 

(ii). In a layered material with a square-symmetric lat- 
tice, the spin-orbit-assisted e-ph coupling is also able to 



generate triplet pairing with definite phases. For sim- 
plicity (and with the reference to an almost cylindrical 
7-sheet in the Fermi surface of Sr 2 Ru04 Q]) we consider 
below a metal with a cylindrical Fermi surface. Sim- 
ilarly to the BCS theory, the effective attraction U in 
Eqs. i|4!7[l is point-like, so that only electrons from the 
same layer can form a pair. Thus, we shall consider only 
two-dimensional order parameters (the interlayer tun- 
neling can be treated using the Lawrence-Doniach ap- 
proach ^|). In this case, triplet p-wave pairs can be 
characterized by the component J z of the total angular 
momentum in the direction perpendicular to the layers, 
which is the sum of the projection L z = ±1 of their 
orbital angular momentum and S z = 0, ±1 of the pair 
spin. For J z = ±1 and J z — ±2, the quantum num- 
ber J z fully determines the order parameter, whereas for 
J z = 0, pairs should be additionally classified accord- 
ing to whether they are symmetric (s) or antisymmetric 
(a) with respect to in-plane x y reflection, for which 
we reserve notations J z — s and J z — a , respectively. 
The second line in Table ITT1 summarizes tensor structure, 
B ms = P^e^D of the order parameter for phases with 
J z = s ,0 a ,±l,±2, where matrices P^ s form the orthog- 
onal basis which realise irreducible representations of J z , 
and where D 2 = J2 ms \B ms \ 2 - 

Using matrices P Jz the effective interaction U can be 
written down in the separable form, Eq. Q) with 



rpSiS-2 _ 



m\m2 J z 



{x^gi+~>cjJl)n^n^Pl^,P 



J, * 



(a c /uj c f p 2 F /2m c , C=-L, ||, 



which is determined by couplings @ to optical phonons 
polarized across (( =_L) and parallel (£ = ||) to the 
layer. The values of coefficients xj z and xj a are listed 
in Table [HJ which shows that positive coupling con- 
stants (determining the attraction in the correspond- 
ing Cooper channel) may appear only for electron pairs 
in the following channels: (a) J z — a where go a = 
K ^gio + ^o«5so = 2g^ - 9so, and (b) J z = ±1 where 
9±i = x±\g~t + &±i0so = fllo - gio- This indicates that, 
within the proposed mechanism, only two triplet phases 
are possible. 

To determine which phase is realized we need to com- 
pare the relative strength of the electron coupling to the 

2 JSC 



in- and out-of-plane sublattice vibrations. If g!! < ^g. 



so that 5o a > g±i, the ground state with J z — 0" is re- 
alized. This phase is characterized by the vector order 
parameter d and an isotropic gap, 



dna = »- 




p - y c p + 

A = 2uj ± e- 2 ' v \ 

g = 2sd - gh- 



(15) 
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TABLE II: Triplet pairing order parameter for electrons with cylindrical Fermi surface in a layered material with square lattice 
and the interaciton constants in the corresponding J z channels due to the spin-orbit-assisted e-ph coupling. 

J z = ±2 J z = ±1 J z = s J z = a 

(ABM) 

P^Z = ±(S mx ±iS my )(5 sx ±i8 sy ) Pt^ = -^ 2 {5 mx ±i8 mv )8 sz Pf m = ^{5 mx b~ sx +5 my 5 sy ) P°^ = ^ (S my 5 SX - S mx 5 ay ) 
x±2 = x±i = — 1 >fo s = —2 HQa — 2 

ii±2 = — 1 ><±i = 1 >co B — — 1 Ao a = — 1 



Here [2(j, D — 2 3 / 2 (lu±/ g)e 2 / vg was found from mini- mizing £ in Eq. iTPTjl with respect to D. 



If 



<7s'o > ffiw so that 5±i > 3o a ( m terms of 
microscopic parameters, for a stronger coupling to 
the in-plane sublattice displacement, / (y/rn^u)^) > 

aj_/(^/TOj_cjj_)), then the ground state belongs to the 
phase with J z = ±1 which is similar to the Anderson- 
Brinkmann-Morel (ABM) phase discussed in the theory 
0,0] of superfluid He 3 . This phase also has a circularly 
isotropic gap: 



i±i 




E p = y /e»+A*, 
A = 2co|| e - 2 /^, (16) 
.9 = gL - gio- 



In conclusion, we have shown that the spin-orbit- 
assisted e-ph coupling can generate triplet pairing with 
a definite tensor structure of the p-wave order parame- 
ter in each type of a crystal. For example, in a layered 
metal where conducting planes have square lattice struc- 
ture, the proposed mechanism of triplet pairing leads to 
two possible phases in Eqs. l|15l) and (|16|l . one of which 
(ABM) is actually considered [l| as the most plausible 
candidate for the triplet order parameter in supercon- 
diting Sr2RuOi that has the layered perovskite struc- 
ture with body-centred tetragonal group symmetry. On 
the one hand, the latter coincidence may be purely ac- 
cidental, so that a detailed analysis of the phonon spec- 
trum and the spin-orbit-assisted e-ph coupling strength 
in each particular material based upon first-principle cal- 
culations is needed. On the other hand, the proposed 
pairing mechanism suggests a n isotope effect of the usual 
sign in the triplet superconductivity materials (for exam- 
ple, using oxigen isotopes in SrzRuOi, as reported in Ref. 
|2ll l22j ). which would manifest the involvement of e-ph 
coupling in the formation of triplet pairing. 
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